Self-localization of composite spin-lattice polarons 
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Self-localization of holes in the Holstein t-J model is studied in the adiabatic limit using exact 
diagonalization and the retraceable path approximation. It is shown that the critical electron- 
phonon coupling Ac decreases with increasing J and that this behavior is determined mainly by the 
incoherent rather than by the coherent motion of the hole. The obtained spin correlation functions 
in the localized region can be understood within a percolation picture where antiferromagnetic order 
can persist up to a substantial hole doping. These results restrict the possibility of self-localization 
of holes in lightly doped cuprates. 
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The interplay of strong electronic correlations and 
electron-phonon (EP) interactions in the formation of 
dressed quasi-particles is one of central puzzles of high- 
Tc superconductors. Recent ARPES experiments in un- 
doped cuprates were interpreted in terms of strong EP 
coupling giving rise to self-localization of holes Q, 13 - 
Similar effects can be expected in the manganites in the 
colossal magnetoresistance regime, where polarons are 
not only dressed by spin and lattice excitations but also 
polarize the orbitals in the neighborhood of the holes. 
The properties of composite polarons are subtle and not 
well understood, as the different mechanisms involved 
may support or compete with each other 0, 0, IS S El ■ 

Electronic motion in weakly doped Mott insulators is 
determined by the constraint of no double occupancy of 
sites which renders the motion of holes predominantly in- 
coherent with an energy scale t, the bare hopping energy 
of electrons. In addition, the antiferromagnetic (AF) ex- 
change interaction J allows for spin flips leading to co- 
herent hole motion at the bottom of the band with a dis- 
persion characterized by J. In the presence of strong EP 
coupling the spin polaron transforms into a spin-lattice 
polaron. The formation of this composite polaron af- 
fects in general both the coherent translational as well 
as the internal incoherent motion. In recent studies of 
Holstein t-J models appropriate for cuprates, it has been 
found that the effect of EP-interaction on spin-polarons 
is strongly enhanced as compared to polarons in uncorre- 
lated systems. In particular, the critical EP-coupling Ac 
for self-localization of composite polarons is significantly 
reduced [11111,11,11,0,111, and in the regime A < Ac the 
coherent polaron mass is strongly enhanced as compared 
to the uncorrelated case 



In spite of several studies there are important open 
questions, which we address in the following: (i) Is the 
critical EP coupling Ac for self-localization determined by 
the coherent bandwidth or rather by the incoherent hole 
motion? (ii) What is the dependence of Ac on J? (iii) 
Furthermore, is the observed low doping concentration 



X ^ 0.02 for the destruction of AF long-range order in 
cuprates compatible with self-localized polarons? 

To answer the above questions we investigate the Hol- 
stein t-J model in the adiabatic limit where the kinetic 
energy of the lattice can be neglected: 
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The first two terms in Eq.l^ represent the Hamiltonian 
Ht-j of the t-J model. The transfer matrix elements 
tij include both nearest and next-nearest neighbor hop- 
pings t and t' , respectively. The third term describes 
the lattice potential H^p proportional to the EP cou- 
pling constant g, the local displacement fields Ui, and 



the hole density operator = 1 



In high-Tc cuprates this term reflects the breathing mo- 
tion of oxygen ions around the hole and the subsequent 
change of the Zhang- Rice singlet energy due to the local 
lattice distortion. The last term is the elastic energy Hph 
with force constant K. In the adiabatic limit displace- 
ments Ui can be treated classically and are in the equili- 
birum determined by Ui = [g/K] < >. At vanishing 
and weak EP-coupling H exhibits coherent (delocalized) 
quasi-particles with bandwidth ~ J. Nevertheless, we 
will show in the following that many features of H can be 
very well represented by the simpler Holstein t-J^ model 
in the retraceable path approximation (rpa), where the 
motion of carriers is entirely incoherent. It is convenient 
to express results in terms of the dimensionless EP cou- 
pling parameter A — g^/8Kt. For noninteracting elec- 
trons we have then Ac ~ 1. In the following we also put 
t = I< = 1. 

The Hamiltonian Eq.(^ is solved by exact diagonal- 
ization (ED) of small planar systems with a square lat- 
tice. The determination of the critical EP coupling is 
straightforward, namely, it separates ground states with 
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homogeneous and inhomogeneous lattice dispacements Ui 
corresponding to a delocalized and a self-localized po- 
laron solution, respectively. Since the itinerant ground 
state for a hole in the t-J model is at ko ~ (7r/2,7r/2), 
the ground state is degenerate. This degeneracy leads in 
general to inhomogeneous Ui which spoil the interpreta- 
tion 1^ 1^ . We avoid this problem by choosing twisted 
boundary conditions which lift the degeneracy. For fixed 
Ui we find then the ground state using ED for systems 
with TV = 18 and 20 sites, whereby the equilibrium is 
reached by the iteration of the self-consistency relation 
= (g/K) <n'^>. 

The t-Jz Holstein model is solved within rpa, which 
excludes loop motion and thus neglects a very small co- 
herent bandwidth that arises at small Jz from loop tra- 
jectories Using a Bethe lattice the expectation value 
of the electronic energy, namely, {Ht-j^) + {Hep), is given 
by the lowest energy eigenvalue cq of the symmetric ma- 
trix Hill with the elements 
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Here I — 0,1,2.. denotes the shell of 1-th nearest neigh- 
bors consisting of Ni ions with Nq = 1 and iV/ = 4 • 3'"^ 
for I ~ 1,2, .., and ui is the displacement of one of the 
equivalent ions in the shell /. The expectation value of 
the density at the ion i, (n'l), is equal to \e{0 , / Ni , 
where e(0, /) is the normalized eigenvector belonging to 
the lowest eigenvalue eg and / is the shell index of the ion 
i. Finally, the total energy cq -I- (Hph) is minimized with 
respect to the displacements {ui}. 

Fig. ^a) shows the electronic energy E^i — {Ht-j) + 
{Hep) and the total energy Etot = Eei + {Hph) as a func- 
tion of A for the Holstein t—J model and J = 0.1, 0.3 and 
0.5. The results were obtained by ED for one single hole 
and clusters with A'' = 20 sites. For A ^ 1 the hole is 
completely localized yielding Eei — 2Etot ~ —8tX. In the 
limit A ^ the well-known single-hole energies E^i of the 
t — J model are reproduced. For J ^ the energy is close 
to the rpa result E^^ ^ —2^/3t (Note that we avoid in the 
ED studies very small values for J where the ground state 
is Nagaoka-type ferromagnetic with E^^ ~ — 4t). For a 
finite J Eel increases because the itinerant hole weakens 
the antiferromagnetic bonds. Using ED results for spin 
correlation functions this increase is about 4.72J which 
agrees well with Figma). For A < Ac the homogenous 
solution with E^^ is stable (the small slope appearing 
in Fig^a) is a finite size effect since Ui ~ 1/TV). For 
A > Ac the localized solution has the lowest total energy 
approaching the linear behavior at large A. Due to the 
first-order transition between itinerant and localized so- 
lutions the ED data exhibit a small jump in Eei and a 
change in slope in Etot at Ac. 




FIG. 1: (color online) Electronic and total energies Eei and 
Etot, respectively, as a function of A for (a) the Holstein t-J 
model using ED and (b) the Holstein t-Jz model using rpa. 
The inset explaining the curves in (a) also holds for (b). 

Fig. n^b) shows the same quantities as in Fig. Q^a) but 
calculated for the Holstein t-Jz model using the rpa. The 
curves are nearly identical with those of Fig. ^a). The 
main difference is the absence of a well defined transition, 
at least in the cases J = 0.3 and 0.5. For J = 0.1 Etot 
and, to a lesser degree, Eei are practically constant up 
to A '--^ 0.4. While for small J Etot decreases gently with 
increasing A Eei shows a sudden decrease at A ~ 0.5 
before reaching its asymptotic linear behavior. Similar 
features are seen in the curve for J = 0.1 in Fig. QJa) 
obtained by ED. In the limit J Etot and Eei are 
identical and equal to — up to a critical coupling 
aI^'' ~ 0.580. For A > A^^'' there exist two solutions of 
the extremal equation. One is the homogenous solution 
where Etot and Eei are identical and independent of A. 
The second solution describes a localized polaron where 
Etot exhibits an upward jump by about 0.4 and Eei a 
downward jump by about 0.6 at A^"^. This means that 
the localized solution is unstable in the interval [Ai"-* , Ac] , 
where Etot of the localized solution crosses Etot of the 
homogenous solution at Ac ~ 0.662. With increasing J 
the jumps are replaced by cross-overs with decreasing 
absolute changes until for J > 0.06 there exists only one 
solution of the extremal equation describing a crossover 
from an extended to a localized polaron. 

Fig. 13 shows ED data (circles and squares) for Ac as 
a function of J. The dashed curve represents data for 
t' = —0.3, i.e., a case with a larger and more dispersive 
coherent part than for t' = 0. As a result, the dashed 
curve lies above the ED data for t' = 0, showing, that 
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FIG. 2: (color online) Dependence of Ac on J comparing 
results from ED, rpa, and Born approximation. 

decreasing the mass M for the coherent motion leads to 
an increase in Ac as expected. Increasing J also decreases 
M but Fig. 121 implies that in this case Ac decreases and 
does not increase. This is unexpected, since an interpre- 
tation in terms of an opposite trend \c J has been 
given in a study of the Holstein t-J model 0. 

In order to understand the dependence of Ac on J we 
have estimated Ac using the rpa and the Born approxima- 
tion. For J < 0.06 the rpa leads to a first-order transition 
between a localized and delocalized ground state, the re- 
sulting Ac's are shown as crosses in Fig. |21 In order to 
find a sharp transition at larger values of J the energy of 
the homogeneous solution has to be compared with the 
localized one. An estimate for the loss of coherent kinetic 
energy in the localized state is AEkin = e — e(ko) ~ 0.65J 
|lO| where e(k) is the quasiparticle dispersion and e its 
average value. Ac then follows from the rpa result via 
the relation AEkin = -E'tot(O) — Etoti^c)- The result- 
ing values (crosses for J > 0.1 in Fig. give the right 
trend but cannot reproduce quantitatively the strong de- 
crease of Ac with J of the ED results. Modifying the Born 
approximation to take inhomogenous local potentials 
into account while retaining the homogenous self-energy 
we have determined the smallest value of A where a lo- 
calized solution exists and identified this value with Ac. 
The resulting diamonds in Fig. El are rather near to the 
ED data. 

The above results may be interpreted in simple physi- 
cal terms as follows: At J = the hole motion is entirely 
incoherent and determined by the energy scale t. The 
rpa solution shows in this region a sharp transition at 
Ac from a very localized polaron to a very extended ob- 
ject. The obtained critical value Ac ~ 0.6 is somewhat 
reduced compared to the free case ~ 0.85. For larger J 
the transition in rpa evolves into a pronounced crossover 
in £'fot(A) as is evident from Fig. 2b), where the crossover 
value A* moves to lower A's due to the reduced polaron 
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FIG. 3: (color online) (a) Density no in the center of the po- 
laron as a function of A. Results obtained by ED for the t-J 
model (symbols) are compared with rpa results for the t-Jz 
model (lines), (b) Log-linear plot of displacement fields ui of 
the Holstein t-Jz model versus distance I showing the expo- 
nential decay of hole distribution in the self-localized regime. 

radius. Although the sharp transition finally involves 
AEkin oc J, Ac is effectively governed by the incoherent 
solution and its A* . The decrease of Ac with J reflects a 
smooth transition from the large energy scale t charac- 
teristic for the incoherent motion to the scale J relevant 
for the coherent part. This picture may also explain ED 
data showing that Ac actually increases with J at large 
J ^ t where the incoherent part is small compared to 
the coherent one. Fig. Ola) shows the hole density at the 
center of the localized state, denoted by riQ, as a function 
of A for different values of J. The symbols are results 
from the ED of clusters with = 20 sites in the lo- 
calized region, the lines correspond to the rpa. no = 1 
corresponds to a completely localized and no — l/N to 
an extended state. The agreement between both results 
in Fig. Oa) is excellent, showing, that the degree of lo- 
calization of the polaron is to a very good approximation 
determined by the incoherent part of the hole motion de- 
scribed well by the rpa. In the limit J — )■ the rpa yields 
a curve for uq which is zero for A < A^'^'', jumps to about 
0.65 at aI"', and then smoothly approaches the J = 0.1 
curve. The crossover value A* is clearly visible in Fig. 
|3Ja) and can be identified as the inflection point in the 
no(A) curves. Fig. ^h) illustrates the qualitative change 
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FIG. 4: Averaged spin correlation function < Sf Sj > for the 
Holstein t-J model using ED, 1 hole, and 18 sites, = 1,4 
denote first- and third-nearest neighbors. 

of the displacements ui at Ac, i.e., for A > Ac decays 
exponentially with distance I. 

An important consequence for magnetism follows from 
self-localization of holes in weakly doped high-Tc super- 
conductors. It is well known that AF long-range order 
vanishes in hole-doped cuprates at quite low doping con- 
centration XAF ^ 0.02 — 0.04. This small value is due 
to the incoherent motion of holes which destroys the AF 
order; on the other hand a percolative model based on 
static holes and broken AF bonds leads to much larger 
critical concentrations xaf ^ 0.5. Self-localization at 
strong EP coupling A > Ac would hinder the incoherent 
motion of holes and preserve the antiferromagnetic order. 



for the considered doping x ~ 0.05 long-range AF order 
persists throughout the localized region and that immo- 
bile holes reduce mainly the magnitude of the SCF's but 
not their spatial decay. 

In conclusion, we have shown that the critical EP cou- 
pling Ac for self- localization in the Holstein t-J model 
decreases with J and that the crossover between basi- 
cally incoherent (coherent) motion of the hole at small 
(large) J is responsible for this decrease. Exact diagonal- 
ization results for Ac, the density and displacement distri- 
bution as well as spin correlation functions have been ob- 
tained and successfully interpreted in the localized regime 
within the retraceable path approximation. Our results 
suggest that for dopings x ^ 0.05 self-localized holes can- 
not suppress sufficiently antiferromagnetic correlations 
to explain the observed absence of long-range order in 
the cuprates at this doping level. We therefore believe 
that holes in lightly doped cuprates are not self-localized, 
and that their EP coupling is below the critical value 
Ac ^ 0.25 for J = 0.3. Our conclusions are consistent 
with the observation of a sharp quasi-particle peak in 
ARPES experiments of 3 % doped LaSrCu04[l3, and 
also agree with the explanation of transport and optical 
data in lightly doped La2Cu04 in terms of holes bound 
to impurities at low but delocalized at high temperatures 

Useful discussions with O. Gunnarsson are acknowl- 
edged. 



The strong enhancement of AF correlations in the self- 
localized regime A > Ac is displayed in Fig. 2]which shows 
ED data for the first {(P — 1) and third {cP — 4) near- 
est neighbor spin correlation function (SCF) of the t-J 
model. For the undoped Heisenberg antiferromagnet the 
corresponding SCF values are -0.116 and 0.067, corre- 
spondingly. For A >> Ac the hole is totally localized and 
one expects a decrease of these values due to the four 
broken bonds yielding a reduction factor of 16/18 for 18 
sites. The resulting values -0.103 and 0.059 agree quite 
well with the numerical data at large A in Fig. ^ sug- 
gesting the applicability of a simple percolation picture 
at large A. The figure also illustrates that for a given 
type of neighbors the SCF become independent of J well 
above Ac which reflects the fact that the localizing EP 
coupling and not the string potential ~ J, Eq.Q, deter- 
mines the extension of the polaron. The slight decrease 
of SCF's with decreasing A in the localized region can be 
understood in the rpa. With decreasing A the hole ex- 
plores more and more neighboring sites producing hereby 
an increasing number of broken bonds which reduce the 
SCF. Such a picture reproduces for A > 0.6 quantita- 
tively the ED data for SCF in the Ising, and to a good 
degree, also in the Heisenberg case. This suggests that 
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